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Abstract 



We explore the applicability of the exact renormalization group to the study of 
tunnelling phenomena. We investigate quantum-mechanical systems whose energy 
eigenstates are affected significantly by tunnelling through a barrier in the poten- 
tial. Within the approximation of the derivative expansion, we find that the exact 
renormalization group predicts the correct qualitative behaviour for the lowest en- 
ergy eigenvalues. However, quantitative accuracy is achieved only for potentials with 
small barriers. For large barriers, the use of alternative methods, such as saddle-point 
expansions, can provide quantitative accuracy. 
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Introduction: The exact renormalization group |T|, |^] is a powerful method with a wide range of 
applications in various fields. It provides a framework in which it is possible to study non-perturbative 
aspects of physical phenomena. This is made possible by an exact renormalization-group equation 

||, |3j that describes the dependence of generating functionals for the correlation functions of the 
theory on a coarse-graining scale. In particular formulations, such as the effective average action 
[||, ||], an exact equation can be obtained for the flow of the coarse-grained free energy of the system, 
a quantity with intuitive physical interpretation. 

The biggest difficulty one faces in this approach concerns the approximations that must be made 
in order to turn the exact renormalization-group equation, a functional differential equation, into 
evolution equations for quantities such as the effective potential. A widely used approximation method 
employs an expansion of the effective action in the number of derivatives of the fields appearing in it. Its 
validity, at the formal and practical level, has been studied and tested extensively |§-[|9]]. The absence 
of a small expansion parameter makes the precision of a given truncation in the number of derivatives 
hard to estimate, and one often must rely on comparisons with alternative methods. However, one can 
obtain answers for quantities that are difficult to compute, such as critical exponents, amplitudes and 
the critical equation of state for second-order phase transitions, or the bubble- nucleation rate beyond 
the semiclassical level for first-order phase transitions ||. 

In this letter we are interested in the applicability of the approach to tunnelling phenomena. 
We work within the formulation of the effective average action §, g. Recent work 0, |l|] has 
demonstrated how the presence of a coarse-graining scale leads to a consistent quantitative description 
of such phenomena, even beyond the leading semiclassical level. Of related interest is the question of 
the convexity of the effective potential in field theory, which is induced by tunnelling configurations. 
This issue has been addressed within the renormalization-group approach as well 12]. 



The study of the above problems within field theory is often obscured by technical difficulties, 
such as the need to regulate the ultraviolet divergences, or possible infrared problems arising from 
the presence of massless modes in the phase with spontaneous symmetry breaking. A much simpler 
framework is provided by quantum mechanics. One does not have to deal with ultraviolet problems 
and spontaneous symmetry breaking does not arise. Moreover, the exact answers are well known 
and comparisons are straightforward. A recent study []l3| has also explored the issues we address in 
this letter. However, it employed an expansion in powers of the field, an approximation that has a 
significant effect on the results. Our treatment goes beyond this approximation and leads to a more 
complete conceptual understanding. 

The problem: We would like to understand how well quantum-mechanical tunnelling is described 
by the exact renormalization-group approach. This approach employs field-theoretical language. How- 
ever, a theory of a real scale field in one time and zero space dimensions can be viewed as a quantum- 



mechanical system. This mapping is obvious in the functional integral formulation [14| through the 
replacement of the field <j) by the position x of the quantum-mechanical particle. As a result, ap- 
proximations employed in the study of field theories in more than one space-time dimensions can be 
checked through comparison with exact results in one dimension. 

We employ the formalism of the effective average action r& ]|, [| that can be interpreted as a 
coarse-grained free energy. All fluctuations of the system with characteristic momenta larger than an 
infrared cutoff k (and wavelengths less than 2iT/k) are integrated out and incorporated in the effective 
couplings appearing in IV The dependence of on the coarse-graining scale k is described by an exact 
renormalization-group equation [|j . An expansion in derivatives of (ft results in evolution equations for 
the functions multiplying the various terms in IV These form an infinite system of coupled partial 
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differential equations in <ft and k. In practice one must truncate this system, by keeping only a finite 
number of terms in the action. 

In one Euclidean space-time dimension, the lowest truncation level (often referred to as the local 
potential approximation) uses an action of the formF] 

(1) 

At the next level the derivative term is multiplied by a non-trival wavefunction renormalization Z^cf)), 
while higher levels include more derivatives of the field. For the approximation of eq. (jl|), the exact- 
renormalization group equation for results in the evolution equation |], 13] 




dk 2ir \ k 2 



for Uk(4>)- This equation has been derived for a sharp infrared cutoff k that acts like a ^-function, 
completely integrating out fluctuations of (j) with momenta q > k and excluding fluctuations with 
q < k. It also appears in other formulations of the exact renormalization group in the sharp-cutoff 
limit |, 

In ref. jO^] Uk(4>) is further expanded in powers of <fi, and eq. (§) is turned into a system of ordinary 
differential equations for the coefficients of the expansion. This infinite system is truncated at a finite 
level. In ref. H it was shown that such a truncation does not reproduce correctly the solution of eq. 
(0) in the regime that is relevant for the convexity of the effective potential. For this reason, we do not 
make this additional approximation here, but solve numerically the full partial differential equation 
(Q) instead. For this we employ numerical algorithms discussed in ref. |15| j. 

The boundary conditions for the solution of eq. @ are fixed at a scale k = A much larger than the 
physical scales of the low-energy theory. At this "microscopic" level, the potential is determined by the 
fundamental theory. The case of one dimension is particular in this respect, because the contribution 
of the ultraviolet regime to the evolution of the potential is negligible. This reflects the absence of 
ultraviolet divergences in quantum mechanics. As a result, one can interpret the initial condition 
Ua as the potential at the level where no fluctuations of the system are taken into account, i.e. the 
classical potential V. The solution of eq. (|2|) in the limit k — > incorporates the effect of fluctuations 
at all scales, and U$ becomes the effective potential U. In quantum-mechanical terms, U determines 
the expectation value of the energy of the system for a given expectation value of the position. 

We consider ^-symmetric classical potentials of the form 

U A = V = -\m^ + \^ + ^. (3) 

Their minima satisfy U(<j) m i n ) = 0. By expressing all dimensionful quantities in units of m, we can set 
m = 1 in the above relation. The presence of a barrier implies that tunnelling configurations can be 
relevant for this system. For A —* the height of the barrier becomes much larger than the distance 
between the minima and the origin, and we expect tunnelling to play an important role. 



The solution: Starting with the above initial condition, we solve the partial differential equation (||) 
numerically. A typical solution is presented in fig. 1 for A = 0.05. The evolution starts at k = A S> 1 
with Ua = V . It is apparent that for k S> 1 the potential Uk changes very little. This means that the 

1 We use h = 1 throughout the paper. 
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ultraviolet (short-distance) fluctuations do not contribute to the evolution, and reflects the absence 
of ultraviolet divergences in quantum mechanics. For k = 0(1) the evolution becomes fast, while 
the location of the minima moves to zero. Eventually, the evolution slows down for k — » and the 
potential takes its final form with only one minimum at the origin. 
Several properties of this solution have physical significance: 

• In one dimension, the location of the minimum always runs to zero for k — > 0. This reflects the 
absence of spontaneous symmetry breaking in quantum mechanics. 

• The value of the potential at the minimum changes from zero to a positive value. In more than 
one dimensions, the absolute scale of the potential is related to the cosmological constant and 
is very sensitive to the contribution from the ultraviolet fluctuations. In one dimension, the 
ultraviolet contributions are negligible and the absolute scale of the potential is meaningful. In 
fact, the value of Uq = U at the minimum corresponds to the expectation value of the system in 
the vacuum, i.e. the eignevalue Eq of the ground state. In our case Eq = U(0). 

• In field-theoretical terms, the first excitation of the system above the vacuum corresponds to a 
particle at rest. Its mass term is given by the second derivative of the potential with respect 
to (f) at the minimum. This means that the eigenvalue of the first excited quantum-mechanical 
state can be computed as E\ = Eq + \/d 2 U (0)/dcp 2 . 

• If the lowest energy eigenvalue Eq is below the maximum of the barrier of the classical potential 
V(0) = 1/(16A) (as in our example), tunnelling plays an important role in the problem. The 
solutions are then of trully non-perturbative nature. 

In fig. 2 we plot [d 2 Uk((j))/d(j) 2 ] jk 2 during the whole evolution. The continuous lines depict the 
initial part of the running with A > k <; 1, while the dashed ones the final part with 1 i> k > 0. We 
observe that there is a range of k for which d 2 Uk{4>) / d(f) 2 ~ —k 2 near the origin of the potential. The 
reason for this behaviour is the presence of a pole at —1 for [d 2 Uk((ft)/d(j) 2 ] jk 2 in eq. (||). As this pole 
cannot be crossed, d 2 Uk{(f))/d(f) 2 remains close to —k 2 while k is reduced towards zero. The pole is 
approached only for A ^ 0.1 and the part of the evolution the system spends near it grows for A — > 0. 
In this limit the importance of the barrier increases. 

This behaviour has been observed in studies of field theories with spontaneous symmetry breaking 
in higher dimensions ]12| . It results in the presence of a flat part of the potential between the two 
minima for k — ► 0. (d 2 Uk^o((f)) / d(j) 2 becomes zero between the minima.) It cannot be reproduced 
through the expansion of Uk(4>) m powers of <f> employed in ref. In the one-dimensional case 

we are studying, the minimum of the potential always moves to the origin and d 2 Uk{ty / d4> 2 becomes 
positive at the end of the evolution. However, the longer the part of the evolution the system spends 
near the pole, the smaller the final positive value of d 2 Uk^o(0) / d(f> 2 . As this value determines the 
difference between the first two energy eigenvalues, the presence of the pole in the evolution equation 
predicts a vanishing E\ — Eq for A — > 0. This is in qualitative agreement with the expectation from 
quantum mechanics. 

The quantitative accuracy of the predictions can be checked through comparison with the results 
of the numerical solution of the Schrodinger equation with a potential given by eq. (^). In fig. 3 we 
plot the values of Eq, E\ predicted by the solution of the evolution equation, along with the exact 
results, as a function of A. We also show the height of the barrier V(0) = 1/(16A). Comparison with 
Eq, Ei determines how important tunnelling is for the corresponding solutions. In fig. 4 we give 
values for the difference AE = E\ — Eq. We also include the prediction of the dilute- gas instanton 
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approximation [jO], [l3j 

For A ^ 0.15 the agreement between the values of Eq and E\ resulting from our analysis and the 
exact results is at the few-per-cent level. For A = 0(1) it becomes better than 1%. The deviations are 
in the same direction, so that our prediction for the difference AE agrees at a level better than 1% 
with the exact value for all A ^ 0.15. We conclude that the lowest order of the derivative expansion 
(eq. (||)) leads to very good quantitative results when the potential does not have a large barrier. This 
conclusion is reaffirmed through the study of a potential without a barrier. In fig. 5 we give the results 
for the unharmonic oscillator, corresponding to the potential of eq. (^) with a mass term m 2 = — 1. 
All this is in very good agreement with ref. [JD1], which indicates that the additional approximation 
of the potential by a polynomial, employed there, is very good for large A. 

For A 0.15 the lowest energy eigenvalue is below the top of the barrier of the classical potential. 
This means that tunnelling plays a significant role in the solution. In the region 0.09 A 0.15 
our results have an accuracy better than 10%. We conclude that the renormalization group accounts 
reasonably well for tunnelling effects in this region. For A ^ 0.09 the energy eigenvalue of the first 
excited state is below the top of the barrier as well. Our results have large deviations from the exact 
values and the renormalization-group approach fails to give a good quantitative picture. However, 
the correct qualitative behaviour is still reproduced. In fig. 4 we observe that the difference AE 
diminishes with decreasing A. Excessive requirements in computer time, in order to reproduce correctly 
the solution near the pole, forbid the numerical integration of eq. @ for very small A. Already for 
A = 0.04 the numerical error becomes close to 1%. However, as we explain below, we expect that 
AE becomes for A — > 0. An extrapolation of our results in this limit, as well as an approximate 
analytical solution (see below), verifies our expectation. However, the correct quantitative dependence 
of AE on A, given by the dilute-gas instanton approximation (eq. (||)), is not reproduced. 

Discussion: In order to understand better the nature of the solutions of the partial differential 
equation (||), it is useful to consider some approximations that permit an analytical treatment. For 
large k the potential has a minimum 4>o(k) away from the origin. The /c-dependence of the minimum 
can be obtained by considering the total ^-derivative of the condition dUk/d^l^^ = 0. One finds 

#o = 1 dU' k (4> ) = 1 U'j!'(M 1 

dk U'^o) dk 2tt U'^o) k 2 

where primes denote derivatives with respect to eft. For small A the potential Uk(4>) can be approxi- 
mated by a quartic polynomial, as in eq. @. Moreover, the coefficient of the quartic term has a weak 
dependence on k, as can be verified through the numerical solution. Within perturbation theory, one 
can explain the smallness of the omitted corrections by the fact that they involve powers of A. 

We consider first the case in which U^((po)/k 2 <C 1 during the whole evolution. Within our 
appoximations, eq. (||) can be integrated easily with the result 

*§(*) = $(A) -~i (6) 

where we have assumed A 3> 1. For any initial (j)o(A) there is a value k cr = 3/(7t</>q(A)) = 12A/7T at 
which the minimum moves to the origin. This is in agreement with our expectation that spontaneous 



1 + 



k 2 



(5) 
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symmetry breaking cannot appear in one dimension. At scales below k cr the minimum remains at the 
origin, while the mass term obeys the equation 

du£(p)_ iu^(o)\ 1+ u^y\-\ 



dk 2tt k 2 



k 2 



For positive uj^(0), the mass term grows as k is further reduced. The evolution stops for very small 
k, when U£(Q) » k 2 (decoupling). 

It is possible that t/^Oo)/*; 2 = 0(1) before the minimum of the potential moves to zero. In the 
approximation of a quartic potential with a constant quartic coefficient A & 0.1, we have U^(4>o)/k 2 ~ 1 
for k p ~ 1 > k cr . This complicates the analytical study of eq. (|5|). Moreover, it raises the possibility 
that, as k 2 is reduced below k 2 , it may become much smaller than U£((j)o) and the evolution may stop 
with the system in the phase with symmetry breaking. This scenario is never realized because of the 
presence of a pole at d 2 Uk/d(j) 2 = —k 2 in the evolution equation (|2|). This pole is approached first at 
the origin, where eq. ([?]) applies. The relevant behaviour was analysed in detail in ref. [12| for field 



theories in more than one dimensions, where symmetry breaking can occur. It was shown that the 
solution approaches the limit k — > with U£(0) ~ — k 2 . The curvature at the origin vanishes for k = 
and this leads to the convexity of the effective potential. Polynomial approximations of the potential, 
such as the one employed in ref. jl^], cannot reproduce this behaviour and lead to the appearance of 
singularities at non-zero values of k. 

Our numerical solution, depicted in figs. 1 and 2, reproduces the correct behaviour for A = 0.05. 
The inner part of the potential becomes very flat as the pole is approached. At some stage the flatness 
approaches the minima and prevents the decoupling behaviour that would occur for U'^(<po) S> k 2 . 
Eventually the minimum moves to the origin and the system settles in the symmetric phase. Because 
of the induced flatness near the origin, the mass term d 2 Uk->o(4>) / d(j) 2 that determines the energy 
difference E\ — Eq is very small. As this type of behaviour becomes more pronounced in the limit 
A — > 0, we expect that E\ — Eq becomes zero in this limit. 

The unsatisfactory aspect of our study is that the correct quantitative behaviour for A — > 
is not reproduced by the renormalization-group approach. This is apparent in fig. 4, where the 
energy difference E\ — Eq is depicted. This quantity approaches zero linearly with A, instead of 



exp 



-1 



/ 



as predicted by the dilute-gas instanton approximation. We conclude that the 
first order of the derivative expansion cannot account quantitatively for the non-local effects associated 
with tunnelling. 

However, the use of an alternative method, such as an expansion around a saddle point, can remedy 
the situation. A very good example is provided by the study of tunnelling rates in three dimensions 



[ 1 1 1 . The renormalization group can be used for the integration of high frequency fluctuations (that 
may even induce radiative symmetry breaking), while an expansion around a saddle point (such as 
the scalar instanton or bounce) can account for tunnelling effects. The complementarity of the two 
methods guarantees the consistency of the calculation. 

An approximate analytical solution: As a final confirmation of the previous discussion we 
present an approximate analytical solution of the evolution equation near the origin of the poten- 
tial. By defining the quantities p = (j) 2 /2, p = kp, Uk(p) = Uk(p)/k, t = log k, we can write eq. (§) 

as 

dt ~ 2u * P U * 2vrl + u>+2~pu>r (8) 



6 



Primes on Uk{p) denote derivatives with respect to p. Near the origin of the potential {p ~ 0) we can 
approximate the above equation by 

du' b I 3 1 \ du' h , , , 

-w + y + ^iT^)ii +2u ^ - (9) 

This is a first-order partial differential equation for Uj.(p), which can be solved with the method of 
characteristics. 

The most general solution is 

P\f U ^k+ 2 -tan_1 VXfc = F ( U 'k e2t ) for U 'k > 

3 




u' k + — log ( | = G ^ e ^ for u / < . (1Q) 

The functions -F and G are determined by the initial condition for the potential at k = A. For the 
choice of eq. (pi) with m 2 = 1 we find 



1 v , — 3 , / 1 + 



G(x) = + log| iili 

The above solution displays all the characteristic behaviour we discussed earlier. The minimum 
of the potential always moves to the origin at a scale k cr = 12\/tt, while the pole at k~ 2 d 2 Uk/d(j) 2 ~ 
u'u = —1 is never crossed. We can also obtain an approximate expression for the mass term at the 
origin (p = 0) in the limit k -> 0, u' k (0) -► oo. We find y/dU k ^o(0)/dp = ^d 2 U k ^ (0)/d(/) 2 ~ 6A, for 
small A. This explains the linear dependence of AE on A, observed in fig. 4 for small A, even though 
the predicted slope is not quantitatively correct. 
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Figures 



• Fig. 1: The evolution of the potential as the coarse-graining scale is lowered from k = A ^> 1 to 
k = 0. The initial form of the potential is given by eq. (]3|) with m 2 = 1 and A = 0.05. 

• Fig. 2: Same as in fig. 1 for [d 2 Uk/dcp 2 ] /k 2 . 

• Fig. 3: The first two energy eigenvalues Eq, E\, as predicted by the exact renormalization group 
in the lowest order of the derivative expansion, along with the exact results, as a function of A, 
for the potential of eq. (|3|) with m 2 = 1. The dotted line indicates the height of the barrier for 
the potential of eq. (||). 

• Fig. 4: The energy gap E\—Eq as a function of A. We display the exact values, the predictions of 
the exact renormalization group in the lowest order of the derivative expansion, and the results 
of the dilute-gas instanton approximation. 

• Fig. 5: Same as in fig. 3 for the potential of eq. (H|) with m 2 = —1. 
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